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We demonstrate that the amplitude {,fn\d" (q"jv'y5q)\Q) does not vanish in the limit of zero quark 
masses. This represents a new kind of violation of the classical equation of motion for the axial 
current and should be interpreted as the axial anomaly for bound states. The anomaly emerges in 
spite of the fact that the one loop integrals are ultraviolet-finite as guaranteed by the finite-size of 
bound-state wave functions. As a result, the amplitude behaves like ~ in the limit of a large 
momentum p of the current. This is to be compared with the amplitude (ttI^" (971^75?) |0) which 
remains finite in the limit — > 00. 

The observed effect requires the modification of the classical equation of motion of the axial- 
vector current by non-local operators. The non-local axial anomaly is a general phenomenon which 
- , is effective for axial-vector currents interacting with spin-1 bound states. 

^ ■ PACS numbers: 11.30.-j, 11.40.-q, 11.40.Ha 
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The analysis of two-photon decays of pseudoscalar mesons in the late 40-s ||] led to the discovery of the famous 
g3 \ axial anomaly the divergence of the axial vector current violates the classical equation of motion and does not 
' vanish in the limit of zero fermion masses. For a quark of mass rriq and charge eQq the properly modified equation of 
motion contains a local anomalous term and has the forni|^ 
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d''{qiul^q) = 2zm,g759 + N.^^^^FF, F^, = e^^^pF"". (1) 
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that the form factor G'*' defined by the 2-photon matrix element 



, with F^i, the electromagnetic field tensor. This modification of the classical equation of motion accounts for the fact 

in 

(7(91)7(92)19'^ (97.759)10) =e2e,,.J,,e*G^(p^g^ 92'), (2) 
QQ ' does not vanish for uiq = but turns out to be a constant independent of the current momentum p = qi ^ q2- 



G^{p\ ql = ql = 0) = -2N,{QqflATT\ (3) 



' In this letter we study the properties of axial currents when one of the photons, 7(92), is replaced by a vector meson 
(~| 1^(92), e. g. a /9-meson. We demonstrate that the form factor defined according to the relation 

(7(9i)l^(92)|a^(9-'7.759)|0)=ee,,,.,,,5G^(p2,9?,9i) (4) 

has also an anomalous behavior and does not vanish for massless quarks. This occurs in spite of the fact that the 
(~| , vector meson is a bound 99 state and the corresponding loop graph has no ultraviolet divergence. Moreover, the 
I!" • anomalous behavior is observed for both, the neutral and the charged axial-vector currents. 
The classical equation of motion reads 



5_j . d^iq'-fy-f^q) = i(mg + m'q)q'-f^q + e{Qq> - Qq)A''q'-ful5q, (5) 

where is the electromagnetic field. The (7T^|...|0) matrix element of the second term on the r.h.s. of vanishes 
to order e. Therefore, to this order, the classical equation of motion (|^) predicts G^ = for niqi — niq — 0. We find 
however to order e and for large 

~ Mvfv/p\ (6) 

where My and fv denote the mass and the decay constant of the vector meson. Because of the dependence on p^, 
the newly found deviation from the classical equation of motion for bound states corresponds to a non-local anomaly. 

We are interested in the region \p'^\ ^ m^,AQ(^^, in which case the quarks in the triangle diagram have high 
momenta and their propagation can be treated perturbatively. We discuss in parallel the 77 and 'yV final states in 
order to show how the anomaly emerges in both cases. As in Ref. we consider the spectral representation for the 
axial- vector current itself before forming the divergence. In distinction to |^ , where the spectral representation in p^ 
was considered, we use the spectral representation in the variable 9I. This allows us to take bound state properties 
into account. 



1 We use the following notations: e = V47raom, 7^ = i7''7S^7^ = -1, Sp (7^7^7''7''7'') = 4ie'"^"'', 

dy.Ai, — d^Afj,. eabcd — eai3iii>a°'y^ c^d" for any 4-vectors a,b,c,d. 
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A. The absorptive part of the triangle amplitude 



The amplitude of the single-flavor axial current between the vacuum and the two-photon and the photon-vector 
meson states, respectively, can be written in the form 



^*''(g2)e*"(gi)r.„M(Zi,'72). 



(7) 



The absorptive part ti^ap of T^afs is calculated by setting the two quarks attached to the external particle with the 
momentum (72 on the mass shell, see Fig 1. 





FIG. 1. Diagrams describing (7(gi)7(g2)lg7i^75(?iO) and/or {'y{qi)V {q2)\q"f,^j5q\0) , with p = qi + q2. The cut corresponds to 
the calculation of the absorptive part in the variable g|. 

t^ap is our basis for the spectral representation of T^ap in terms of the variable (/f. The coupling at the vertex f3 is 
ipe-Qq if the particle 2 is a photon, and — 7/35(92 )/'\/^ if i* is ^ vector meson.^ The couphng will be further 

discussed below. By taking the trace and performing the integration over the internal momentum in the loop it is 
straighforward to obtain ti,cep- The result is automatically gauge- invariant 

9?tiya/3(gi,g2) = 0, q2tual3{qi,q2) = 0. (8) 

It is therefore possible to write the covariant decomposition of t^a,p{qi, (72) in terms of three invariant amplitudes 

)ZC2. (9) 

This Lorentz structure is chosen in such a way that no kinematical singularities appear. We take 7(91) to be a real 
photon, qf — 0. Hence, the term containing the invariant amplitude ci does not contribute to the divergence of the 
current. Setting in addition niq — one obtains for cq and C2 with s = 



47r (s — p2 



C2(p^s) 



An 



(10) 



where ({s) = 2NcQ'l9{s) for the 77 process and C,{s) 



iq^x"/ / ; J--"--™'- ^v"y — -'2\/N2Qqg{s)9{s) for the process. 

Clearly, the absorptive part t^afsiqi^ q2) of the axial- vector current matrix element respects the classical equation 
of motion, that is 



p^co(p^,s) - SC2{p'^,s) = 0. 



(11) 



B. The triangle amplitude and its divergence 

The full amplitude T^afsiqi, 92) has the same Lorentz structure as its absorptive part 

+ {ql ^vfiaqi q2p^uq2Ciqi )iC2. (12) 

The absence of any contact terms in T^afi can be verified by reducing out one of the photons and using the conservation 
of the electromagnetic current. The invariant amplitudes Ci can be represented by the following dispersion integrals 



^ The full vertex has the form Q F^ = —l/^/Nc [jis — ^^2m (^^ ~ fe);3/(-\/s -I- 2m)]g{{ki + fe)'^), but the term proportional 
to (fci — k2)i3 does not contribute to the trace. The overall (— ) sign is the standard choice of the phase of the vector meson 
wave function which leads to a positive leptonic decay constant. 
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C.,{p\ql^Q^ql) = - r^^lfLds. (13) 
Jo s- q^-iQ 



All the integrals converge and thus need no subtraction. 
Taking the divergence of T^af3 we find 



tp T^ap = <P / —ds-q2 —ds\ tq^aq^fS- (14) 

L "'o * ~ 92 "'0 ^ 92 



The form factor G defined in Eqs. (^ and (^) now reads 

In the case of the 77 process C,{s) is a constant. The integral can be performed and gives the well-known value shown 
in Eq (|^). In the case of the matrix element the integrals converge even better since g{s) which appears in C(s) 
descibes the spatial size of the vector meson. We conclude from Eq. (|l5|) that the divergence of the axial-vector 
current is nonzero for triq ^ not only for the 77 but also for the final state! Namely, 

G^'ip^ql) - 2VA^eQ,-4 T^^ds. (16) 



The behavior with respect to p^ is however different from the 77 case and has the form (p^) ^ l/p^ for the large 
values of p^ where our formula applies. 

For the transition to the 7p (isospin-1) arising from the isovector axial current we obtain 

The parameter k in this equation is non zero for rriq = and jp^j 00. fp is defined by the relation {p~\d'y^u\0) = 
fpMptl. 

Eq. ( p^ ) takes into account the soft contribution to the form factor Qp . For large jp^j one should take care of the 
QCD evolution of the p- meson wave function from the soft scale /i^ ^ 1 GcV^ to the scale /i^ ~ This can be 

done most directly by expressing k in terms of the p-meson light-cone distribution amplitudes |^ 

1 1 

(p((Z2)M»7A^(0)|0) = ~iq2x{e*x)fpMp j due'^'i-^Hu) + eypMp j due''-'''^g'^l\u), 



1 

(p((Z2)M»7A75^(0)|0> = -^ex,,p^e*^qPx''fpMpJ due'^''''^ g^l\u) . (18) 



The diagrams of Fig 2 explain the procedure we follow: The quark propagator connects the axial current with the 
photon and the two distant space time points are bridged by the p meson wave function. This leads to the following 
expression for k 

y 1 

■qi 





u 



FIG. 2. The diagrams for the matrix element (7(gi)p^(g2)|d7iy75f |0), with two space time points bridged by the p meson 
light-cone distribution amplitude. 



^We want to point here to the similarity of the form factor with the tt'j transition form factor F„jj» (p^)- For a detailed 
analysis of the latter we refer to Ref. Likewise, the form factor C'' describing the amplitude (ppld" {q~fv'y5<l)\Ci) has some 
common feature with the pion elastic form factor. 
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du 







1 + U (a) I - U (j,) 



-9T' 1 — * 



(19) 



The leading-twist distribution amplitudes g'f ^ = 6m(1 — u), g^""* = | (l + (2u — 1)^), and <f> = — m)(2m — 1) give 
the main contribution for large and lead to the value k ~ —3/2. Corrections to this value are calculable in terms 
of the higher twist distribution amplitudes 

Summing up, our results are as follows: 
1. The divergence of the axial vector current and thus the form factor G does not vanish in the limit — > 0. This 
holds for the 7p final state as well as for the 77 final state. The observed effect for vector mesons requires a proper 
modification of the equation of motion for the axial current. For large momenta p of the axial- vector current the 
corresponding 'bound state anomaly' can be described in terms of a non-local operator appearing at order e (there 
are no local operators of the appropriate dimension): 



Single flavor current : 
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^''(97.759) = 2zm q^^q + If^i^F/' + en Q, {d^{q^,q) ■ F^-^} + 0{a-^). (20) 
Isovector charged current : 

du[d'-i^-fzu) = i(m„ + md)u-fzd - ie((3„ - Qd)u'y^j5d ■ A" + eK-^^^^-^-^^D"^ ^df_,{uj^d) ■ F^^j + 0{a^^). (21) 

K can be expanded in a power series of as- In leading order one obtains the value k = —3/2. 

As pointed out in [Q, the anomaly is important in rare decays: for instance, in i? — > p7 decays, it substantially 
corrects the weak annihilation amplitude which carries the CP violating phase. 

2. The amplitude for the pp final state also stays finite for mq = 0. The corresponding non-local anomalous term for 
the divergence of the isovector axial current d,j{qj^j5T"'q) appears already at order 0(e^). The operator structure 
of the anomalous term is more complicated in this case. One of the possible lowest-dimension operators which has a 
nonvanishing (pp\...\0) matrix element and thus will contribute to the anomalous term (in leading order in l/p^) is 
the product of the two isovector tensor currents 

^a6CgM-/3n-l (qa^,T\ ■ qao^pT-q) . (22) 

Accordingly, for large \p'^\ the {pp\...\Q) amplitude of the divergence of the axial- vector current is given by the factorized 
matrix element of the anomalous term and has a 1 /p^ suppression 

{p{qi)p{q2W{d^.^,um = e,,..,,.^^^ GPP ^ fl/p" + 0(l/p^). (23) 

3. We have illustrated the appearance of non-local anomaly due to vector mesons in QCD. This anomaly is of a 
general nature and should be present in any theory containing = 1^ bound states. For example, the anomaly 
will also contribute to the generation of orthopositronium by the leptonic axial- vector current. In contrast to the 
conventional local anomalies, there is no obvious cancellation of the non-local anomalies in the standard model. 
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